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Abstract 

The density-matrix and Heisenberg formulations of quantum mechanics follow — for unitary 
evolution — directly from the Schròdinger equation. Nevertheless, the symmetries of the corre- 
sponding evolution operator, the Liouvillian C = i [•,-£/], need not be limited to those of the 
Hamiltonian H. This is due to C only involving eigenenergy differences, which can be degenerate 
even if the energies themselves are not. Remar kably, this possibility has rarely been mentioned 
in the literature, and never pursued more gener ally. We consider an example involving meso- 
scopic Josephson devices, but the analysis only assumes familiarity with bàsic quantum mechanics. 
Subsequently, such C-symmetríes are shown to occur more widely, in particular also in classical 
mechanics. The symmetry's relevance to dissipative systems and quantum-information processing 
is briefly discussed. 
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I. INTRODUCTION 



Density matrices streamline quantum statistics, by combining the calculation of quantum- 
mechanical expectation vàlues and classical averaging over a probability distribution of 
states, into a single trace operation [ï], |2|. The aesthetically inclined might further ap- 
preciate that for a pure state \ip), its density matrix \ip)(ip\ is uniquely defined, unlike \ip) 
itself [|]. 

While the state vector obeys the Schròdinger equation id t \ip) = H\ip) {h= 1), the density 
matrix evolves according to the von Neumann (or quantum Liouville) equation 

d t p = i[p,H)=£(p), (1) 

where £ is the Liouville (super-)operator. Because £, is defined in terms of the Hamilto- 
nian H, the two operators are closely related; in particular, any symmetry of H is inherited 
by C (see Section [TT] for a precise statement). The density-matrix formulation is also uniquely 
useful in more general cases where the Liouvillian no longer is given by C = i [ • , H] @, |5| , or 
where the evolution even is effectively nonlocal in time ||. In such situations, a pure state 
p = p 2 will in general evolve into a mixture, so that a Schròdinger formulation is not possible. 
This occurs when the system being studied is coupled to a "bath," of which complete knowl- 
edge is neither possible nor desirable. Tracing out the bath degrees of freedom then yields 
the reduced density matrix of the system. The latter's effective evolution can deviate from 
unitarity by (a) dissipation (with the attendant fluctuations fTj): changes in the diagonal 
(in the energy basis) elements of p representing occupation probabilities, and (b) dephasing: 
random fluctuations of the system's energy levels, which suppress off-diagonal elements of p 
even without bath-induced dissipative transitions between those levels. Studying and con- 
trolling the non-unitary evolution of p is for instance crucial if its off-diagonal elements are 
used to store (quantum) information ||. 

For the moment, however, we restrict ourselves to unitary evolution (^j. Then, the 
Schròdinger and von Neumann equations would seem to be equivalent, a complete set of 
solutions for the latter being given by 

\tpj)(ijJk\ exp{i(u k -Uj)t} , (2) 

where H\ipj) = It is our purpose to highlight the fact that the reverse need not 

hold, and to examine the consequences. Namely, consider a system such that 

Co>3 — Lü\ = UI4 — LÜ2 7^ . (3) 

For any a, b, the superposition ) (^2 1 +^1 ^3) ("^4 1 ] exp{i(u;2 — ^i)t} is then a solution of ([[]), 
which is not of the form (0) if ah 7^ 0. In other words, the concept of degeneracy is wider for 
the von Neumann than for the Schròdinger equation, since only energy differences have to 
be equal §| — physically appealing, as it is only those differences which are gauge invariant 
and, hence, observable. The generalization to more than four states will be obvious. 

II. FORMALISM OF LIOUVILLE SYMMETRY 

Before discussing specific instances, let us outline the theory of the above-mentioned 
Liouville symmetry or C-symmetry, in a form general enough to apply to classical mechanics 



We will not dwell on the mathematical subtleties which arise in the case of continuous spectra. 
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as well. First of all, a few notions from operator theory should be recapitulated, and some 
notation established. However, readers familiar with superoperators, or on the contrary, 
interested more in applications than in formalism, can skip to Section fTTT] and consult the 
present section as needed. 

Quantum mechanics deals with operators on a Hilbert space of states EL These operators 
can themselves be seen as points in a Liouville space L [|TÜJ, which can be given the inner 
product 

(A, B) = Tr{A^B} . (4) 

As always, the trace can be calculated in any basis, though it has the form Ylji^jl ' iV'j) 
only in orthonormal ones. Sidestepping the issue that the trace may well diverge in infmite 
dimensions, the inner-product axioms are readily verified. In particular, there are no prob- 
lems with density matrices, which are positive-semidefinite operators normalized according 
to (p, 1) = 1. Observables are Hermitian operators A, so that (A) = (p,A) G R. 

Suppose the subspace V Ç L is \-closed: V = VL Then W = {p G V : p = p^} is a 
real subspace spanning V; in particular, observables have real products (Q). An orthonormal 
basis of W then is a Hermitian (C-)basis for V; cf. (^0]) below. To proceed, consider VEI Ç EI. 
Counterintuitively, this need not itself be a subspace; cf. 2 
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which is not in the range of any matrix in the span of the two indicated ones. If however 
VEI in fact is a subspace of EI, one may define its projector P = Pvh- If furthermore 
one has V = PhP, we call V a conventional Liouville subspace, generated by the Hilbert 
subspace VEI. In the notation developed below, V = P\P r for the projector onto V. In 
all other cases, V will be said to be an unconventional Liouville subspace. Note that the 
dimension of a conventional Liouville subspace is a square, and that a proper Liouville 



subspace containing an invertible element is always unconventional; compare both with ( |2Ü|) 
in Section |V|. 

The logical next step is to consider linear superoperators^ S in Liouville space [O]. Or- 
dinary operators A can be promoted to superoperators in several ways. Presently we only 
introduce A\{p) = Ap and A T (p) = pA 



C 



lÜf , so that (JI]) amounts to 



(6) 



A superoperatorS is a left (right) multiplication iff\/A,B : S(A)B = S(AB) [BS(A) = 
S{BA)]. In this case, S = S(l)\ [S = S(t) T ]. For the proof of the reverse implication, 
note that S(t)\(A) = S(t)A = S(tA). The other parts are analogous and/or trivial. 
Elementary calculation rules are [A\, B\\ = [A, B\ and [A T , B T ] = [B, A] v , while [Ai, B T ] = 0; 
also, (A\/ r y = {A^)y x so that one can write A\, without confusion, and li = l r = 1. 

In this language, the quantum-mechanical superposition principle and probability con- 
servation imply that time evolution should have the form p(t) = U(t,t')[p(t')} for t > t! . 
This leads to the theory of completely positive trace-preserving superoperators [IT2L [13 



2 Foregoing Hermiticity, there are also 3D counterexamples. 

3 This term here merely means "operators acting on operators," and has nothing to do with conversion 
between bosons and fermions. Well, almost nothing — see p. O. 
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In its full scope, this paper advocates studying the symmetries of U in Liouville 
space, in any physics problem. 

Explicitly, symmetries are superoperators S commuting with U, 

[U, S} = 0. (7) 

Of course, the commutant of U can always be written down formally in terms of the latter's 
Jordan normal form [T2j . The point is to first find symmetries on physical grounds, and then 
to use these in the spectral analysis of IÀ instead of relying on brute force. This conceptually 
paral·lels the case of conventional Hamiltonian symmetries, though technical aspects may 
well differ. 

Specializing to systems without memory (at most Ohmic damping) leads to a differential 
formulation 4 p = £{t) (p), and the symmetries can be studied on the level of the generator 5 £, 
hence the name Liouville symmetries. If further the system is closed then £ is given by (Q, 
regardless of whether £ and H are time dependent or not. The evolution in Hilbert space 
is now unitary, and £ is readily checked to be anti-Hermitian, 

(A,£(B)) = -(£(A),B) . (8) 

When time evolution is generated by a Hamiltonian, the question arises how the symme- 
tries of £ relate to those of H: operators 6 A such that [H, A] = 0. One has [£, S] = for 
a superoperator S of the form S = A\ or S = A T , and an £ of the form (Qj, iff [H, A] = 0. 
In other words, an ordinary operator can be an £-symmetry only if it is an ordinary sym- 
metry in the first place. In fact, by @ and the calculation rules, [£, A\/ T ] = i[A,H]y x 
is immediate. One can try to extend this result to more general symmetries S. How- 
ever, while [H\, S) = [H r ,S] = =>- [C,S] = by (H), the reverse does not hold (v.i.). 
Clearly, £-symmetry is not conveniently studied in terms of the Hamiltonian, even for 
closed systems — commutation with both H\ and H T is too strong, while commutation with 
only one is too weak because e.g. [H\, A r ] = does not imply any symmetry. 

To give some crèdit to H, its very existence will in general lead to simple £-symmetries. 
Namely, in terms of the eigenstates H\%pj) = U)j\il>j), we can first of all define the projec- 
tor Vjk(p) = PjpPk (i-e., Vjk = Pj,iPk,r), where Pi = \ipi){ipi\ are ordinary Hilbert-space 
eigenprojections. 7 Note that Vjk is an £-symmetry and commutes with both H\ and H r 
(immediate from the calculation rules), but is not an ordinary operator. Slightly more in- 
terestingly, one has the transfer %k{p) = \ipk)(ipj\p\ipj)(4>k\, satisfying [£,Tjk] = — any 
mixture of eigenstates is stationary. On the other hand, [Hy r , Tj^\ = iff ujj = u^. By defi- 
nition, it further follows that £-symmetry carries one ordinary symmetry to another. That 
is, if [H, A] = then also [H, S(A)] = 0. Of course, there is no guarantee that anything 
nontrivial ensues for any given system, i.e., S(A) could well be in the span of {1, H, A}. 



4 In e.g. quantum òptics, one often speaks of a master equation. 

5 Also at this level, there is a correspondence between Schròdinger \(A) — (C(p),A)] and Heisenberg 
[(A)' = (p,tf(A))] pictures. 

6 Note that A need not be Hermitian, cf. angular-momentum raising/lowering operators for spherical sys- 
tems, and most unitary operators. 

7 Projectors corresponding to degenerate eigenspaces of H are readily written down as a sum of such Vjk- 
For instance, V H = ^ w ._ Vjk performs an energy measurement. 
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We call a superoperator S real if it maps Hermitian p to Hermitian p'. From this it follows 
that S(p*) = [S(p)y for general p. The evolution operators LoxlÀ are necessarily real for any 
physical system (in fact, stronger conditions hold); on the other hand, iC is only Hermitian 
for unitary evolution. In contrast, A\ is Hermitian for any observable A by the calculation 
rules, but real only if A = al, a G M; apparently, the two are independent concepts. Note 
now that, e.g., the theorem relating C- and ií-symmetries would not have come out so neat 
if one had considered S = A Y A\ instead of S = Ay T , cf. the examples of Section [VT|. Thus, 
we have not demanded our £-symmetries to be real. This seems reasonable, given that 
eigenvectors (|2|) of C are non-Hermitian if they correspond to a nonzero eigenvalue, cf. the 
complex-conjugate eigenmodes of the underdamped classical harmònic oscillator. That is, 
one cannot formulate the whole theory in terms of Hermitian p only. 

The above abstraction from "degenerate energy differences" to "commutation with £" is 
good for, if nothing else, arriving at the classical limit. Classical mechanics is defined on 
an even-dimensional phase-space manifold T |Ï4] , on which observables / are real functions. 
Technically, it is optimal to start from a symplectic 2-form. For our purpose, however, 
it is sufficient to have a Poisson bracket 8 f,g i— > {/,<?}, and local canonical coordinates 
(p, q) = (pi,..., p n , qi,..., q n ) in terms of which 

dqj dpj dpj dqj 

The Poisson bracket is then verified to satisfy the Jacobi identity, and one can further define 
the phase-space volume dpdq, which is invariant 9 under any canonical transformation, i.e., 
(p, q) >— > (p ; , q') preserving the form of (||). Define the inner product (/, g) = f r f*g dpdq, and 
a density p as a positive phase-space function. Then, normalization of p and the calculation 
of expectation vàlues proceed exactly as in the quantum case. 

By analogy with the preceding, superoperators are defined as operators in the function 
space of observables. Given an observable /, the obvious associated superoperator reads 
g i— > fg. Because of the commutativity of classical mechanics, there is no need to notationally 
distinguish the two roles of /. A superoperator S corresponds to an observable iff S(f)g = 
S(fg), in which case S = SI. The proof is obvious by comparing to the quantum case. 

Another important class of superoperators is derived from mappings A : T — > T, which we 
assume to be invertible. Then we can define 10 g i— > Ag = g o A^ 1 . Such a map is said to be 
canonical if the associated coordinate map is a canonical transformation. More elegant is the 
equivalent invariant criterion V/, g : {Af,Ag} = A{f,g}. A superoperator S corresponds to 
a phase-space map iff the homomorphic property S(f)S(g) = S(fg) holds, and moreover S 
is invertible. Note that the map being canonical is a separate condition. To prové this, one 
can reconstruct the mapping by taking sequences of functions approaching a 5-function on 
phase space, and study their behavior under S; we will not dwell on the technical details. 11 



8 There will be no confusion between the Poisson bracket and the anticommutator, since the former (latter) 
only arises in classical (quantum) mechanics. 

9 Even the orientation of phase-space volume is preserved. 

Note a difference with the quantum case, where the "position" part of such a mapping would induce an 
ordinary operator tp(q) ip(A~ 1 (q)). For instance, parity is a quantum but not a classical observable. 

1 Indeed, the whole phase space can be reconstructed from the observable àlgebra |Ï5|] , so that a superop- 
erator preserving this àlgebra has to preserve phase space. 
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Hamiltonian dynamics is generated by an observable H, according to d t p = {H, p} = Cp. 
An observable / is said to be an integral of H if {H, /} = 0. A canonical map A is said to 
be a Hamiltonian symmetry if AH = ií. 12,13 As in the quantum case, a superoperator S is 
said to be an £-symmetry if [£, S] = 0. An observable is an C- symmetry iff it is an integral 
of H; a canonical map generates an C-symmetry iff it is a Hamiltonian symmetry. For a 
proof, it suffices to write out [£, f] = {H, /} and [C, A]p = {H—AH,Ap} respectively [in 
the former case, (|9]) shows that Poisson brackets satisfy a product rule]. Note furthermore 
that the commutator [H, S] simply does not occur in the classical theory. 

As an afterthought to the above, we point at the possibility that superoperators could be 
declared points in a new space themselves, with hyperoperators mapping between them, and 
so on at infinitum. At least in the finite-dimensional case, the construction should proceed 
without difficulty. It is not clear to us to what extent this hierarchy has been pursued in 
mathematics. 
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FIG. 1: A loop consisting of a d-d junction and a SET, threaded by an external flux. 



III. LIOUVILLE-SYMMETRIC QUBIT CIRCUIT 

Let us give an example, where £-symmetry occurs naturally and actually facilitates prac- 
tical calculations. Consider the circuit in Figure [I], containing an intrinsically degenerate 
phase qubit, coupled to a single-electron transistor (SET) [|TB[ acting as a switch fT7| . Read- 
ers not interested in Josephson devices can skip the circuit analysis and simply take its 
effective Hamiltonian ( |T2D below as a postulated one, acting in a 4D Hilbert space. The 
exact nature of the qubit is immaterial here; it can for instance be realized as a d-d grain 
boundary. A loop is formed by coupling the SET to the bulk of the qubit. The SET is taken 
superconducting as well, with Josephson couplings E 1;2 , so that phase coherence can be 
maintained throughout the loop, and a phase frustration e can be imposed by an external 



12 Note that the invariance A{H) = H did not play a role in the quantum case 

13 It is well known that a smooth family of such symmetries generates an integral of H — in effect a finite- 
dimensional Noether charge. However, there does not seem to be a counterpart for £-symmetries, since 
these are linear superoperators already so that taking the infinitesimal generator of a smooth family 
gives nothing new. Compare to quantum Hamiltonian symmetries [ií, ^4(A)] = 0, where passing to the 
generator, [ií, 9aA|a=o] — 0, merely yields another Hamiltonian symmetry. 
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magnètic flux. The SET and qubit phases are S and 0b respectively, with conjugate charges 
Q s and Qh- The SET island can be polarized by a gate voltage V g , coupled through a gate 
capacitor C g . If the loop inductance is negligible and the SET (qubit 14 ) capacitances 
(C*b) satisfy Ci <C C^, Standard electrostatics yields the Hamiltonian 

H = H Q + E h ((j) h ) - Ex cos(0 b -0 s ) - E 2 cos(0 s -0 e ) , (10) 

- 2 ^ Qs+C ^> U/(0,C7 E ) C? ) [ Qs+ C g V g 

where C s = C x + C 2 + C g . 

The qubit Josephson coupling i?b(0b) is not specified in detail, except by stating that 
it has two degenerate minima at 0b = ±0o <C 27r. If phase tunneling between them is 
small but not negligible, it is then appropriate to go to the qubit basis 0b = 0o0"b, so that 
the free qubit Hamiltonian reduces to H h = Ql/(2Cb) + E h ((f> h ) \— > — Acrg, with A the 
tunneling amplitude. The remaining qubit operator occurring in ( |Ï0|) becomes i— > qa^, 
with q = i(— 0o|Qb|0o) ^ ^- Note that q, being determined by tunneling, is of the order 
of A, while e.g. sin0 o can be considerably larger. The derivation of a finite-dimensional 
effective Hamiltonian may be completed by making the two charge-st&te approximation for 
the SET, i.e., Q s i-> Q + e(l+a*), upon which e ±i<t>s h-> af = (of ± i(T%)/2. Elementary 
manipulations lead to 

ClÇ<5 Ei COS 0o+^2 cos e £ 2 sm e eQ 

fíeff = -Acr b + ^T^b g CTs 2 s ~C~ a * 

Cxqe Eisin0o , . 

~CCT h a 2 ^ ^ ^ 

= a<r£ + + 7 < + í<7* + tat + <>K + wlal , (12) 

where Q = Q + e + C g V g is the effective charge offset from degeneracy. 



IV. HAMILTONIAN ANALYSIS 



Focus on the case Q = 0, where the SET is degenerate, i.e., acts as an open switch. 
By Eq. flï2[) , this implies f3 = e = 0. 15 Then, the quartic equation for the eigenenergies uo 
reduces 16 to a quadratic equation for u 2 , with solutions 

u 2 = o? + 7 2 + S 2 + ( 2 + r] 2 ± 2v/(a 7 -Cr/) 2 + 5 2 (a 2 + V 2 ) ; (13) 

labeling < co> 2 < cj 3 < co> 4 , these have precisely the difference property @. The corre- 
sponding eigenvectors [in the basis \cr^at) = (|++), |H — ), | — h), | ))] are 

/ UUJ n \ 

«[7 + i{f]+5)} + v n (iÇ-a) 
u(i(-a) + v n [y + i{i]-5)] 



\n) 



\ 



(14) 



14 Cb also contains, and in fact may be dominated by, the qubit-SET lead capacitance. 

15 In fact, the vanishing of e is not essential, so that the analysis could be done for any Q if the bit-SET 
charge coupling oc C\q is negligible. 

16 This happens whenever (3íQ = 0. 
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lüZ. Their norms are 



where u = 2(«7 — (r/ — iaS) and v n = a 2 + 7 2 + (í+r?) 2 + £ 2 
(n|n) = Aulvnivn - 25r]). 

The pairing of states (\n), |n)), with 0% = — u; n , is readily explained by {H, S} = for 
S = 0b<7g. For the vectors flU]), one has S'In) = — \n). 

To study qubit evolution in the presence of the SET, one can define the correlator (a° = 1, 
3 t^O) 

/^(í) = TrK(í)a b Vf} • (15) 
Introducing the parity of the Pauli operators 



one has 



implying 



Sa J h/s = p>oÍ /R S 
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f í+\ TV J q2 iHt j —iHt k £\ 

= p j p k p í Tr {Se- iHt a J h e iHt ay s S} 

= p'pVfM-t) , 



(16) 



(17) 



fjki(t) 



E 

Here, (a) nm = (n\a\m) and u mn = uj m - u n . Easy properties are fju(0) = 46j k 5to, \ f\ < 4 
in general, and /^(O) = if [cr^, a*] = 0. To dispose of trivialities upfront, note that since 
the identity matrix does not evolve, fjoo(t) = Trcr^ = 0. 

Of the remaining 45 cases, an exhaustive search, using computer àlgebra to evaluate (|Ï8|) 
in terms of the states (13)— (H3), shows that no less than 25 cancel: 



E 



(n|n)(m|m) 

Im[(cr^) nm (agq-f) mn ] 
(n|n)(m|m) 



cos(w mn t) , p'pV 



1 ; 



l t) , p>p k p í = — 1 . 



(18) 
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(19) 



While 



can be simplified further using, e.g., {aQ 



[oi 



KOnn 



KOnn 



h)nn V b/nn 

[by the conventional symmetry A = Srja^ + (07+^77) cr^Ug + 
aSa^a^ — a deformed SET charge conjugation], these do not account for ( |T9"D in a system- 
atic, unified way. Apparently, the above expansion in eigenstates, while yielding explicit 
expressions for the nonzero correlators, is not suited to explain the vanishing of the other 
ones. 



V. LIOUVILLIAN ANALYSIS 



The above is to be contrasted with the Liouville approach. Let us return to the full 
Hamiltonian (|Ï2|) , and in line with its representation H = J2jk^jk a b ® °~s a l so write p = 
Ylijk Pjk a h ® a s- This generalizes the widely used Bloch vector for two-level systems flSfl . Of 



course, poo = 1 Tr p is conserved; for the other elements, the Standard commutation relations 
for the Pauli matrices yield (IÏÏ) as 
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The antihermiticity of the operator C implícit in fl20|) is an easy consequence of the com- 
mutation with a Hermitian H in (^), see (]8|). For p — > 0, one sees that C decomposes into 
an upper 5x5 and a lower 10 x 10 block. This immediately explains the cancellations of / 
[and many others, if one e.g. also studies Tr {c^W^ÍO^b^s}] n °ted below (|T8|). In NMR, 
a Liouvillian analysis is sometimes also known as the "direct method," especially when de- 
termining the spectrum of C; the correspondence of the latter's eigenvalues to excitation 
freqüències bears some resemblance to Green's function methods ||. While block-diagonal 
Liouvillians have been observed previously [jTjJ, these have seldom been pursued as due to 
a bona fide physical symmetry, which is not necessarily linked to damping or decoherence. 



Let us momentarily return to the detailed Hamiltonian (|lT|) . For E\ = 0, it describes a 
flux qubit coupled to a fluctuating charge. An advantage of flux qubits is that they can be 
designed insensitive to charge noise. Here, weak decoherence is due to the "transition dipole 
moment" q. For an interpretation, note that the qubit capacitor will have a voltage across 
it during flux jumps, which may well couple to neighboring charges; cf. [p0[| . Presently, 
the usual charge-noise insensitivity argument only predicts that the qubit stays coherent 
if Ei — q — 0, a correct but trivial statement for H as in (|TT|) with only a, 7,5, e ^ 



(see Section [Vï] for more on uncoupled subsystems). We stress that our discussion, leading 
to (p0"D, concerns a different phenomenon; cf. note [15]. 

One experimental consequence may be formulated as follows. Start with the qubit pseu- 
dospin polarized in the y-z plane, p(0) = |(1 + a^rcosO + a^r sin9) <g> p s (0 < r < 1). If at 
any later time t one measures this spin along any direction in the same plane, the outcome 
distribution is independent of the SET's initial state p s , even though the two systems may 
interact strongly and anisotropically. This statement, rigorously proven above, does not 
seem to follow from any conventional consideration of this simple two-spin system. 

In detail, the £-symmetry uncovered above reads 

[C,V] = Q (21) 
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[cf. (0)], where V = (V£ — V®) 2 is the projector onto the subspace spanned by the flrst 
five components of p (plus the trivial normalization poo) i n (p0|)- Explicitly: 17 V£(p) = 
|crg Cg> Trb{o"gp} etc. Inside each block, £ can of course be diagonalized separately; this 
would not have been possible using an eigenbasis of the form (|2|), (Jf4]). 



VI. FURTHER EXAMPLES AND DISCUSSION 

In Sections [D []-|V|, the discussion has focused on density matrices and the von Neumann 



equation. However, (2Ï|) shows that what really matters is having an operator (there a 
projector V) commuting with the Liouvillian, as is also apparent in the general treatment 
of Section O. Therefore, the whole discussion can be equivalently formulated in terms of 
the Heisenberg equation, A = —C(A) for any observable A. This shows that the relevance 
of £-symmetry is not in any way limited to quantum statistics. 

Likewise, attention has been implicitly restricted to time-independent Hamiltonians. Now 
suppose that the parameters in (|Ï2|) have an arbitrary time dependence, except for j3(t) = 0. 
The derivation of (EO) remains vàlid, so that all cancellations flT9| ) still hold. However, the 
simple formulas flT8| ) no longer apply, so that in the Schròdinger approach, even brute-force 
analytical verification of these cancellations is hampered. 

Thinking of this situation as interaction with a classical external field, one is naturally 
led to the next generalization. Namely, if the system is coupled to a quantum-mechanical 
"outside," then in the Markov approximation this is reflected by the presence of nonunitary 
terms in C, so that a Hamiltonian counterpart does not exist from the outset. The important 
point is that such a dissipative Liouvillian d may still preserve the block structure of (pÜ|), 
in which case the £-symmetry will survive even in the presence of damping. Examples 
can be immediately given by replacing some of the parameters (other than j3) in (|Ï2"D with 
e.g. oscillator-bath operators, so that the corresponding qubit-SET operators become error 
generators. It is not relevant for the argument whether the reduced evolution in qubit-SET 
Liouville space fl20|) is Markovian or not. Physically, we think of a situation in which the 
circuit Hamiltonian retains its form (fï0| ) (at least to leading order in some damping pa- 
rameter), but with e.g. e susceptible to (quantum) noise. Intuitively, quantum information 
stored in one of the blocks does not "leak" to other ones, comparable to the behavior of 



quantum error- correcting codes [|2T|. Pending a comprehensive theory of £-symmetric open 
systems, the above re-opens the question whether the notion of decoherence-free subsystems 
indeed is "ultimate" 

Of particular interest would be the case when is nonzero in only one of the blocks, 
so that the density-matrix components in the other block do not feel its influence at all. 
This latter block could then rightly be called a decoherence-free Liouville subspace, in which 
quantum information can be stored error-free, provided that the unitary (infernal) part of 
L is under sufhcient control. In fact, it can be precisely indicated where our discussion 
deviates from, e.g., the — correct — analysis in p3| . The latter's theorems 1 and 2 state, 
in the present terminology, that to every decoherence-free Hilbert subspace Cn one can 
associate a decoherence-free Liouville subspace A4n] we ask if there are M. not of this form. 
Note that the basis vectors of M. may well be "off- diagonal elements" in the usual parlance. 



17 Contrast with the projector describing a measurement of crg, viz., Ai^(p) = h(p+ a^pa^), which retains 
both the crg- and the Éb-components. 
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Therefore, we have insisted on VCd = C&P = 0, where V denotes the projector on .M. In 
other words, it is not sufficient to have 18 VCd = C-dP oc V; this is to be contrasted with the 
non-semisimple case in [ 2J] . First of all, however, it needs careful investigation whether this 
is possible in principle for time evolution which ultimately is Hamiltonian in the "universe" 
of system plus environment. That is, can there exist error generators which vanish in some 
unconventional Liouville subspaces Ai but not in others; cf. (^), where all parameters other 
than f3 occur in both blocks. In summary, one may anticipate interesting future research 
into £-symmetric dissipative models. 

The harmònic oscillator (HO) figures in almost any physics text, and the present is 
no exception. It is described by H — |(p 2 + q 2 ) = a^a + |. In the notation of Sec- 

tion [TT], set S + = a\a r so that S- = a x a\ = S\_. Then [£,5+] = [C, S-] = 0, while 
e.g. [H\,S + ] = S + 7^ 0; all of these are conveniently evaluated in superoperator nota- 
tion, using only the representation (^) of C and the calculation rules, together with the 
Standard [H, a] = —a. S± can also be studied in the context of the coherent M, T£ 



and other interesting non-diagonal states [25| which the HO is known to possess. For in- 
stance, <S_(|a)(a|) = |a| 2 |o;)(a| for the coherent state \a) = exp{cW— a*a}|0). One also 
verifies [<S + ,«S_] = iC, so that the Lie àlgebra of (<S + ,«S_,£) is isomorphic to the one of 
(q,p, 1). Note further that a, aJ are said to be SUSY operators in supersymmetric quantum 



mechanics ||26|| . In fact, the HO with these operators is the prototypical system. We thus 
arrive at a connection between our superoperator symmetry and the conventional "súper" : 
the HO is its own superpartner because deleting the ground state yields the same spectrum, 
apart from a shift by u=l — all its energy differences are degenerate. 

This example is of special interest because the classical limit can be studied in detail. The 
Liouvillian is £ = qd p — pd q , and the system's integrals are of the form / o H for any real /. 
The Hamiltonian symmetries are e^ H ^ c for any real 0, including parity e nC . Observing 
that above S + added one quantum of energy to the system, we set d H = pd p + qd q , upon 
which [£, d H ] = is readily verified. However, obviously [H, d H ] ^ 0. This £-symmetry is 
nontrivial: it expresses that nearby trajectòries, not on the same constant-energy contour, 
keep moving in phase. The transformations generated by d H are not canonical, since they 
do not preserve phase-space volume. Explicitly, A : (p,q) i— > (a(H)p,a(H)q) is associated 
with an £-symmetry for monotonic et. 19 This sheds considerable light on the £-symmetry 
phenomenon: in quantum mechanics, only energy differences correspond to excitations and 
quantum beats, hence are observable. In classical mechanics, only the Hamiltonian flow is 
observable, and £-symmetries preserve it. 

Thus, £-symmetry also plays a role in classical mechanics. Given that in quantum me- 
chanics we contemplated the symmetry's survival under certain types of dissipation, one 
should look for a classical counterpart to that as well. In fact, for harmònic problems 
such a counterpart is well known to be additíve noise and velocity-independent friction, 
which preserve the linearity of the equations of motion. For instance, [C tot ,d H ] = for all 



However, the corresponding condition for the Hilbert-sp&ce projection P and error generators F a , when 
these exist for conventional A4n, will automatically lead to our stricter criterion for Notice further 
that in Liouville space we impose conditions both on PCd and on C^P', in contrast, the Hilbert-space 
generators are Hermitian, so that always PF a = F a P. 

One can consider wider possibilities by not requiring A to be invertible. Allowing a singularity at the 
origin of phase space, exòtics such as a(H) = H^ 1 are also possible; all of these lead to [C, A] = 0. 
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FIG. 2: Phase-space flow for the damped-oscillator Liouvillian £ tot = qd p —pd q +jd H , with 7 = (a) 
and 7 = (b) respectively. Bold arrows indicate the action of the £-symmetry d H , which maps 
integral curves onto each other in both cases. 



C tot = C + •ypdp + fiqd q , 20 which includes both physical Ohmic friction (/i = 0) and the ad 
hoc £ to t = £ + 7<9 H for fi = 7, in which case the phase-space trajectòries become simple 
logarithmic spirals; see Figure |2j. The eigenfrequencies read uj = ±a/1 — (7— /i) 2 /4— -(7+//). 
As a result of scaling symmetry, i.e., linearity, the problem remains tractable considerably 
beyond the case of a single oscillator [ETI, and can be given a description strikingly similar 



to the conservative case |28] . In view of the above, one arrives at the prediction that this sit- 
uation has a parallel for certain ncmharmonic systems, provided that the damped/diffusive 
Liouvillians involved are compatible with the £-symmetries of the conservative dynamics. 

A variation on the HO is the Stark ladder, which has its Hilbert space spanned by 
{\n) : n G Z}. In terms of the shift a = Yl n \ n )( n ^\ satisfying ao) = a)a = 1, the 
Hamiltonian can be given as H = ^2 n n\n)(n\ — A(a + a)). Again [H,a] = —a, so that 
the defmition of S± and their commutators with C, Hi/ r can be immediately transcribed 
from the HO; however, presently [S+, <S_] = 0. The Stark ladder is prototypical of many 
problems with a source, including sòlids in an elèctric field and current- (voltage-)biased 



Josephson junctions in the phase (charge) representation [16]. In particular, if the "shifting" 
degree of freedom is shift-invariantly coupled to other ones, a nontrivial problem ensues. 
While S± then still are exact £-symmetries, these are difficult to exploit using conventional 
methods |2£| . 



The simplest possible case concerns two uncoupled subsystems. Then, one has EI = 
l·lW ®H( 2 ) for the Hilbert space and H = <g> 1 (2) + for the Hamiltonian. The 

spectrum can now be labeled as Uik = to^ + 00^ Vi, j, k, i : U{k — íOu — ^jk — ^>ji [cf. (H)]. 



20 Howcvcr, of the finitc transformations A considered in the previous paragraph, only those with a(H) 
const are permissible if 7 and/or /i are nonzero. 
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It is a staple of physics that finding such a decomposition of (H, H) constitutes a dràstic 
simplification of the problem. However, apparently it only leads to a proper symmetry on 
the Liouville level, as it is energy differences not energies which are degenerate. If S^ 1 ' 2 ^ are 
£-symmetries in the subsystems, then 21 S(p^' (g> p^) = S^(p^) ® S^(p^) is one in the 
composite system. The case of most interest is = 1 and the "forgetful" = Tr. The 
latter can be calculated without any knowledge about the second subsystem, and erases 
all reference to it. While for a single system this would merely be (proportional to) the 
projection on p = 1, presently it leaves a meaningful problem for the first subsystem. We 
shall not elaborate on the classical counterpart, which is broadly similar. 



VII. CONCLUSION 

We have introduced the concept of Liouville symmetry, of which we know no previous 
systematic and general treatment in the literature. Examples are neither excessively rare 
or contrived, nor is the symmetry's actual use unwieldy. Namely, one still deals with linear 
evolution and eigenvalue equations, albeit in a space larger than Hilbert space; for classical 
mechanics, it is Standard to visualize the Liouville flow in phase space. In particular, there 
can be advantages to a Liouville formulation even for closed, unitary systems. We have 
only touched upon the application to open systems and quantum-information processing. 
For instance, the issues of storing, manipulating, and retrieving such information in/from 
the unconventional subspaces concerned, and the extent to which this is at all possible, 
are important and cannot be ignored. Also, group-theoretic aspects await elaboration. We 
hope that the present work may provide both the framework and a stimulus for such more 
detailed investigations. 

In any case, all information stored in the matrix elements of p is physical, and can be 
retrieved through, e.g., quantum-process tomography (QPT) pOj . For the Hamiltonian (|Ï2|). 



this should be eminently feasible using present-day technology, though perhaps sooner for 
nuclear spins than for high-T c Josephson devices. With some urgency, we therefore propose 
an experiment which does exactly that: QPT for an £-symmetric system. 
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